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Solution of Optimal Power Flow Problems using Moment Relaxations 
Augmented with Objective Function Penalization 


Daniel Molzahn, Cedric Josz, Ian Hiskens, and Patrick Panciatici 


Abstract — The optimal power flow (OPF) problem minimizes 
the operating cost of an electric power system. Applications of 
convex relaxation techniques to the non-convex OPF problem 
have been of recent interest, including work using the Lasserre 
hierarchy of “moment” relaxations to globally solve many OPF 
problems. By preprocessing the network model to eliminate low- 
impedance lines, this paper demonstrates the capability of the 
moment relaxations to globally solve large OPF problems that 
minimize active power losses for portions of several European 
power systems. Large problems with more general objective 
functions have thus far been computationally intractable for 
current formulations of the moment relaxations. To overcome 
this limitation, this paper proposes the combination of an 
objective function penalization with the moment relaxations. 
This combination yields feasible points with objective function 
values that are close to the global optimum of several large 
OPF problems. Compared to an existing penalization method, 
the combination of penalization and the moment relaxations 
eliminates the need to specify one of the penalty parameters 
and solves a broader class of problems. 

I. Introduction 

Determining the most efficient operating point for a power 
system requires solution of the optimal power flow (OPF) 
problem. This problem optimizes power system operation 
subject to both network equality constraints and engineering 
limits. Non-linearity of the constraint equations generally 
makes the OPF problem non-convex and can lead to local 
optima [1]. Many solution techniques have been proposed, 
including successive quadratic programs, Lagrangian relax¬ 
ation, heuristic optimization, and interior point methods [2], 
[3]. Some of these techniques calculate at-least-locally opti¬ 
mal solutions to many large problems. However, while many 
local solution techniques often find global solutions [4], they 
may fail to converge or converge to a local optimum [1], [5]. 

Recent research attention has focused on convex relax¬ 
ations of the OPF problem, which provide lower bounds on 
the optimal objective value and can certify infeasibility. Fur¬ 
ther, a convex relaxation based on semidefinite programming 
(SDP) globally solves many OPF problems [6]. 

The SDP relaxation of [6] has been generalized to a family 
of “moment relaxations” using the Lasserre hierarchy [7] for 
polynomial optimization [8]-[10]. The moment relaxations 
take the form of SDPs, and the first-order relaxation in 
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this hierarchy is equivalent to the SDP relaxation of [6]. 
Increasing the relaxation order in this hierarchy enables 
global solution of a broader class of OPF problems. 

The ability to globally solve a broader class of OPF 
problems has a computational cost; the moment relaxations 
quickly become intractable with increasing order. Fortu¬ 
nately, second- and third-order moment relaxations globally 
solve many small problems for which the first-order relax¬ 
ation fails to yield the globally optimal decision variables. 

However, increasing system size results in computational 
challenges even for low-order moment relaxations. The 
second-order relaxation is computationally intractable for 
OPF problems with more than ten buses. Exploiting network 
sparsity enables solution of the first-order relaxation for 
systems with thousands of buses [11], [12] and the second- 
order relaxation for systems with about forty buses [10], [13]. 
Recent work [13] solves larger problems (up to 300 buses) 
by both exploiting sparsity and only applying the computa¬ 
tionally intensive higher-order moment relaxations to specific 
buses in the network. Other recent work improves compu¬ 
tational tractability using a second-order cone programming 
relaxation of the higher-order moment constraints [14]. 

Solving larger problems using moment relaxations is often 
limited by numerical convergence issues rather than com¬ 
putational performance. We present a preprocessing method 
that improves numerical convergence by removing low- 
impedance lines from the network model. Similar methods 
are commonly employed (e.g., PSS/E [15]), but more exten¬ 
sive modifications are needed for adequate convergence due 
to the limited numerical capabilities of current SDP solvers. 

After this preprocessing, the moment relaxations globally 
solve several large OPE problems which minimize active 
power losses for European power systems. Directly using the 
moment relaxations to globally solve more general large OPE 
problems with objective functions that minimize generation 
cost has been computationally intractable thus far. 

To solve these OPE problems, we form moment relax¬ 
ations using a penalized objective function. Previous lit¬ 
erature [16], [17] augments the SDP relaxation [6] with 
penalization terms for the total reactive power generation and 
the apparent power loss of certain lines. Eor many problems, 
this penalization finds feasible points with objective function 
values that are very close to the lower bounds obtained from 
the SDP relaxation. Related work [18] uses a Laplacian- 
based objective function with a constraint on generation cost 
to find feasible points are very near the global optima. This 
paper analyzes the physical and convexity properties of the 
reactive power penalization. 


There are several disadvantages of the penalization method 
of [17]. This penalization often requires choosing multiple 
parameters. (See [18] for a related approach that does not 
require choosing penalty parameters.) Also, there are OFF 
problems that are globally solved by the moment relaxations, 
but no known penalty parameters yield a feasible solution. 

We propose a “moment+penalization” approach that aug¬ 
ments the moment relaxations with a reactive power penalty. 
Typical penalized OFF problems only require higher-order 
moment constraints at a few buses. Thus, for a variety of 
large test cases, augmenting the moment relaxation with 
the proposed single-parameter penalization achieves feasible 
solutions that are at least very near the global optima 
(within at least 1% for a variety of example problems). The 
momentH-penalization approach enables solution of a broader 
class of problems than either method individually. 

This paper is organized as follows. Section HI] introduces 
the OFF problem. Section |III| reviews the moment relax¬ 
ations. Section HVl describes the low-impedance line prepro¬ 
cessing. Section |V] discusses the existing penalization and 
the proposed momentH-penalization approaches. Section |Vl| 
demonstrates the momentH-penalization approach using sev¬ 
eral large test cases, and Section IVIII concludes the paper. 

II. Optimal Fower Flow Froblem 

We first present an OFF formulation in terms of rectangu¬ 
lar voltage coordinates, active and reactive power injections, 
and apparent power line-flow limits. Consider an n-bus 
system, where A/" = n} is the set of buses, Q is 

the set of generator buses, and £ is the set of lines. Let 
+jQDk represent the active and reactive load demand 
and 14 = Vdk + jVqk the voltage phasors at each bus k G 
Af. Superscripts “max” and “min” denote specified upper 
and lower limits. Buses without generators have maximum 
and minimum generation set to zero. Denote the network 
admittance matrix as Y = G -f jB. 

Define a convex quadratic cost of active power generation: 

fck (Vd, Vq) = Ck2 ifpk (Vd, Vq))^ + Cklfpk (Vd, Vq) + CfcQ (1) 

The power flow equations describe the network physics: 

n 

Pok =fpk {Vd, Vq) = Vdk Y, (GikVdi - BikVqi) 

n 

+ 'y ^ (BjfcVdi + GikVqi) + Pok (2a) 

i=l 

n 

Qck ~fQk {Vd, Vq) = Vdk ^ ' ( BikVdi GikVqi) 

n 

+ Vqk'Y^{^ikVdi — BikVqi) + Qok (2b) 

i=l 

Define a function for squared voltage magnitude: 

(14)" = fvk {Vd, Vq) = Vl + Vql (3) 

We use a line model with an ideal transformer that 
has a specified turns ratio : 1 in series with a 11 

circuit with series impedance Rim + jXim (equivalent to an 
admittance of gim+jbim = r——) ^nd shunt admittance 
jbsh.im- The line-flow equations are 


Plm = fplm {Vd, Vq) = (14" + Vq]) ^ 

Pm 

(y^dl^dm SIH (^^m) Qlm COS j'^lm 

(VdzV^gm VqlVdm') (,Qlm Sill (^Zm) COS (^Zm)) j'^lm 

(4a) 


Pml = fpml {Vd, Vq) = {Vdm + Vqrn) 9lm 

(V^zV^TTi V^zV^m) {,Qlm COS (^Zm) ~\~ Sill (^Zm)) j'^lm 
(y^dl^qm (fflm Sill (^Zm) ^Im COS (^Zm)) j'^lm 

(4b) 


Ql^ = fqi^ {Vd, Vq) = - (Fj + Fj) [blm + /tL 

-f {VdlVdm + VqlVqm) {hm COS {9lm) + gim Sin {9lm)) I Tim 

{VdlVqm VqlVdm) {dim COS {9lm.) blm, SlU {9lm.)) jTim, 

(4c) 


Qml = fQml {Vd, Vq) = - {VL + V^m) {blm + 

(VdzV^dm V^zV^m) (^Zm COS (^Zm) Qlm Sill (^Zm)) jTlm 
P ( VdZ^gm COS (^Zm) ^Zm SlU (^Zm)) f'^lm 

(4d) 

{Slm)^ = fsim {Vd, Vq) = {fplm {Vd, Vq)f + {fqim {Vd, Vq))^ 

(4e) 


{Sml)^ = fSml {Vd, Vq) = {fpml {Vd, Vq))^ + {/q^I {Vd, Vq))^ 




(4f) 

The classical OFF problem is then 



min fck{Vd,Vq) subject to 

(5a) 

PGt<fpk{Vd,Vq)<P^r 

MkGM 

(5b) 

QGk<fQk {Vd,Vq)<Q^T 

VkeAf 

(5c) 

{VrY <fvk{Vd,Vq)<{Vf^n^ 

VkeAf 

(5d) 

fsim{Vd,Vq) < {srj:Y 

V {1, m) G £ 

(5e) 

fsml {Vd,Vq)<{ST^Y 

V {1, m) G C 

(5f) 

Vql — 0 


(5g) 


Constraint (5g i sets the reference bus angle to zero. 


III. Moment Relaxations 

Since all constraints and the objective function are poly¬ 
nomials in the voltage components Vd and Vq, the OFF 
problem Q is a polynomial optimization problem that can be 
solved with tools from algebraic geometry. We first review 
the “moment” relaxations from the Lasserre hierarchy [7] 
for polynomial optimization problems and then summarize a 
method for exploiting network sparsity. 


A. Review of Moment Relaxations 

Folynomial optimization problems are a special case of 
“generalized moment problems” [7]. Global solutions to 
generalized moment problems can be approximated using 
moment relaxations that are formulated as SDFs. For poly¬ 
nomial optimization problems with bounded variables, such 
as OFF problems, the approximation approaches the global 
solution(s) as the relaxation order increases [7]. The first- 
order relaxation in the Lasserre hierarchy is equivalent to 
the SDF relaxation of [6]; higher-order moment relaxations 
generalize the SDF relaxation of [6]. 






We begin with several definitions. Let the vector x = 
[Vdi Vd 2 ■ ■ ■ yqn] ^ contain all first-order monomials 
of the decision variables in Q. Given a vector a G N^" 
representing monomial exponents, the expression = 
^di ^d 2 '' ■ defines the monomial associated with x 

and a. A polynomial g (i) is 

9{x)= (6) 

where is the scalar coefficient corresponding to 5:“. 

Define a linear functional Ly {p} which replaces the 
monomials a;“ in a polynomial g{x) with scalars j/q,; 

Ly {9} = SaVa ( 7 ) 

Apply Ly {p} to each element of a matrix argument. 

Consider the vector x = [V^i Vd 2 ^", 2 ] ^ containing the 
voltage components of a two-bus system, where the angle 
reference ( [5^ is used to eliminate Vqi, and the polynomial 
g{x) = - (0.95f + fv 2 {yd, yq) = - (0.95)2 ^ ^ y 2 ^_ 

(Constraining g (x) > 0 forces the voltage magnitude at 
bus 2 to be greater than or equal to 0.95 per unit.) Then 
Ly {9} = - (0.95)2 yoop+yp 2 o+yoo 2 - Thus, Ly { 5 } converts 
a polynomial g (x) to a linear function of y. 

Define the vector x^ composed of all monomials of the 
voltage components up to the specified relaxation order 7 : 

^ [ 1 ydi ... yqn yl ydiyd2 ... 

• yqi yi yiiy<i2 yqiV (s) 

The moment relaxations constrain moment and localizing 
matrices. The symmetric moment matrix M-^ {y} has entries 
ya corresponding to all monomials up to order 27 : 

'^1 {y} - Ly [x^x'l^] (9) 

Symmetric localizing matrices are defined for each con¬ 
straint of ©. The localizing matrices consist of linear 
combinations of the moment matrix entries y. Each poly¬ 
nomial constraint of the form / (x) — a > 0 in © (e.g., 
fv 2 (x) — ^ 2 ™'” ^ 0 ) corresponds to the localizing matrix 

M^-/3 {(/ (x) -a)y} = Ly {(/ (x) - a) x^_px'l^_p} (10) 

where the polynomial / has degree 2/3. Example moment 
and localizing matrices for the second-order relaxation of a 
two-bus system are presented in (fT4l) and (fTSl) . respectively. 

The order -7 moment relaxation of © is 


min Ly < fck ^ subject to (Ha) 

^ I fees J 



Vfc € W 

(lib) 


Vfc G W 

( 11 c) 


Vfc G W 

(lid) 


Vfc G W 

(lie) 


Vfc G W 

(Ilf) 



VkeAf 

dig) 


V (Z, m) G C 

(llh) 


V (Z, m) G £ 

(Hi) 

{y} t 0 


(Hj) 

2/00...0 = 1 


(Ilk) 

t/o.. . 007 ) 0 ...0 = 0 

77 = 1,... ,27 

(HI) 


where ^ 0 indicates that the corresponding matrix is positive 
semidefinite. The constraint (II Ikl) enforces x^ = 1. The 
constraint (II 111) corresponds to the angle reference ( [5^ ; the 
77 in (EB is in the index n + 1, which corresponds to the 
variable Vqi. The angle reference can alternatively be used 
to eliminate all terms corresponding to Vqi to reduce the size 
of the SDP, as in (O-iO. 

The objective function and apparent power line flow con¬ 
straints are quartic polynomials in Vd and Vq. Eor 7 = 1 , use 
a Schur complement formulation of these polynomials [13]. 

The order -7 moment relaxation yields a single global 
solution upon satisfaction of the rank condition 

rank (M-y {?/}) = 1 ( 12 ) 

The global solution x* to the OPE problem © is then 
determined by a spectral decomposition of the diagonal block 
of the moment matrix corresponding to the second-order 
terms. Specifically, let /r be a unit-length eigenvector cor¬ 
responding to the non-zero eigenvalue A from the diagonal 
block of the moment matrix corresponding to the second- 
order monomials (i.e., [M.^ {yY\ 2 -k 2 fc’ where k = 2 n-|-l and 
subscripts indicate entries in MATLAB notation). Then the 
vector V* = \fX {gi-.n + jL{n+i)-. 2 n) is the global solution. 

A solution with rank (M.y {?/}) > 1 indicates that the 
order -7 moment relaxation only yields a lower bound on the 
objective value. Increasing the relaxation order will improve 
the lower bound and may give a global solution. 

B. Exploiting Network Sparsity 

The matrices in the moment relaxations quickly grow 
with both the relaxation order and the system size. Eor 
an n-bus system, the number of rows and columns in the 
order -7 relaxation’s moment matrix is ( 2 n + 7 )!/ (( 2 n)! 7 !). 
Solving second- and higher-order moment relaxations of 
problems with more than about ten buses requires exploiting 
network sparsity [10], [13]. Eull details are excluded from 
this summary for brevity; see [13] for a complete description. 

We use a matrix completion theorem [19] to exploit 
network sparsity. A symmetric matrix W with partial in¬ 
formation (i.e., not all entries of W have known values) 
can be completed to a positive semidefinite matrix (i.e., the 
unknown entries of W can be chosen such that W ^ 0) if 
and only if certain submatrices of W are positive semidef¬ 
inite. Consider a chordal extension of the power system 
network graph0 The matrix W is positive semidefinite if and 

* A chordal extension contains all links in the network as well as additional 
links such that every cycle of length four or more nodes has an edge 
connecting two non-adjacent nodes in the cycle. A chordal extension can be 
calculated using a Cholseky factorization of the network Laplacian mah'ix 
using an approximate minimize-degree permutation to maintain sparsity. 



X2=[l Vdi Vd2 Vq2 Vdi VdiVd2 VdiVq2 Vd2 Vd2Vq2 Vq2 [Note; (|^ is used to remove Vgi] 


(13) 


M 2 {y} = Ly {x 2 xl} = 


Mi{(/v2-(0.95)')t/} = 

J/020 + 2/002 — ( 0 . 95 )^ 1/000 

1/120 + 1/102 — ( 0 . 95 )^ 1/100 
1/030 + 1/012 — ( 0 . 95 )^ 1/010 
1/021 + 1/003 — ( 0 . 95 )^ 1/001 


1/000 

2/100 

2/010 

yooi 

2/200 

1/110 

2/101 

2/020 

2/011 

2/002 


t/100 

2/200 

2/110 

2/101 

2/300 

y2io 

2/201 

2/120 

2/111 

2/102 


1/010 

2/110 

2/020 

2/011 

2/210 

2/120 

2/111 

2/030 

2/021 

2/012 


1/001 

2/101 

2/011 

2/002 

2/201 

1/111 

2/102 

2/021 

2/012 

2/003 


t/200 

2/300 

2/210 

2/201 

2/400 

2/310 

2/301 

2/220 

2/211 

2/202 


2/110 

2/210 

2/120 

2/111 

2/310 

1/220 

2/211 

2/130 

2/121 

2/112 


2/101 

2/201 

2/111 

2/102 

2/301 

2/211 

2/202 

yi21 

2/112 

2/103 


2/020 

2/120 

2/030 

2/021 

1/220 

1/130 

2/121 

1/040 

2/031 

2/022 


2/011 

2/111 

2/021 

2/012 

1/211 

1/121 

1/112 

2/031 

2/022 

2/013 


2/002 

2/102 

2/012 

2/003 

2/202 

2/112 

2/103 

2/022 

2/013 

2/004 . 


2/102 - 

( 0 . 95 ) 

^ 2/100 

2/030 + 2/012 

-( 0 . 95 )" 1/010 

1/021 

CO 

0 

0 

-f 

- ( 0 . 95 ) 

2/202 - 

( 0 . 95 ) 

A 

2/200 

2/130 - 1 - 2/112 

- ( 0 . 95 )" 1/110 

1/121 -f 2/103 

- ( 0 . 95 ) 

2/112 - 

( 0 . 95 ) 

^2/110 

2/040 - 1 - 2/022 

-( 0 . 95 )" 1/020 

2/031 + 2/013 

- ( 0 . 95 ) 

1 

CO 

0 

( 0 . 95 ) 

^2/101 

2/031 - 1 - 2/013 

-( 0 . 95 )" 1/011 

1/022 + 2/004 

- ( 0 . 95 ) 


(14) 


(15) 


only if the submatrices of W corresponding to all maximal 
cliques (i.e., the largest completely connected subgraphs) of 
the chordal extension are positive semidefinite. 

The matrix completion theorem enables decomposition of 
a positive semidefinite constraint for a single large matrix to 
constraints on many smaller matrices. This eliminates many 
terms which do not appear in the constraint equations of the 
OFF problem 0 . See [11]-[13], [20] for details. 

Exploiting network sparsity enables solution of the 
second-order relaxation for problems with up to about 40 
buses. Recognizing that the first-order relaxation is suffi¬ 
cient for large regions of typical OFF problems, previous 
work [13] proposed an iterative algorithm to selectively apply 
the higher-order relaxation constraints to specific regions of 
the network [13]. This algorithm used a “power injection 
mismatch” heuristic to determine where to apply the higher- 
order moment constraints. 

We next summarize this heuristic. A moment relaxation 
yields the global optimum if the moment matrix has rank 
one, in which case the algorithm terminates. If the rank 
of the moment matrix is greater than one, we calculate 
the closest rank-one matrix using an eigen decomposition. 
The power injection mismatch heuristic compares the power 
injections corresponding to the closest rank-one matrix with 
the power injections corresponding to the higher-rank mo¬ 
ment matrix. Typical OFF problems result in small power 
injection mismatches at the majority of buses with a small 
fraction of buses having large mismatches. Each iteration 
of the algorithm in [13] tightens the relaxation by applying 
higher-order relaxation constraints for the two buses with 
greatest power injection mismatch. 

IV. Freprocessing Fow-Impedance Fines 

By exploiting sparsity and selectively applying the higher- 
order constraints, the moment relaxations globally solve 
many OFF problems with up to 300 buses. Solution of larger 
problems with higher-order relaxations is typically limited 
by numerical convergence issues rather than computational 
concerns. This section describes a preprocessing method for 
improving numerical properties of the moment relaxations. 


Fow-impedance lines, which often represent “jumpers” 
between buses in the same physical location, cause numerical 
problems for many algorithms. Fow line impedances result 
in a large range of values in the bus admittance matrix Y, 
which causes numerical problems in the constraint equations. 

To address these numerical problems, many software pack¬ 
ages remove lines with impedances below a threshold. For 
instance, lines with impedance below a parameter thrshz 
are removed prior to applying other algorithms in FSS/E [15]. 

We use a slightly modified version of the low-impedance 
line removal procedure in FSS/E [15] 0 Fow-impedance lines 
are eliminated by grouping buses that are connected by lines 
with impedances below a specified threshold thrshz. Each 
group of buses is replaced by one bus that is connected to all 
lines terminating outside the group. Generators, loads, and 
shunts (including the shunt susceptanes of lines connecting 
buses within a group) are aggregated. The series parameters 
of lines connecting buses within a group are eliminated. 

Removing low-impedance lines typically has a small im¬ 
pact on the solution. To recover an approximate solution to 
the original power system model, assign identical voltage 
phasors to all buses in each group and distribute flows on 
lines connecting buses within a group under the approxima¬ 
tion that all power flows through the low-impedance lines. 

A typical low-impedance line threshold thrshz is 
1 X 10“^ per unit. However, the numerical capabilities of 
SDF solvers are not as mature as other optimization tools. 
Therefore, we require a larger thrshz = 1 x 10“^ per unit 
to obtain adequate convergence of the moment relaxations. 
This larger threshold typically introduces only small errors 
in the results, although non-negligible errors are possible. 

Matpower solutions obtained for the Folish [21] and 
most FEGASE systems [22] were the the same before and 
after low-impedance line preprocessing to within 0.0095 per 
unit voltage magnitude and 0.67° voltage angle difference 
across each line. Operating costs for all test problems were 
the same to within 0.04%. The 2869-bus FEGASE system 

^Lines with non-zero resistances are not considered to be “low 
impedance” by PSS/E. We consider both the resistance and the reactance. 













had larger differences: 0.0287 per unit voltage magnitude 
and 1.37° angle difference. A power flow solution for the 
full network using the solution to the OPF problem after 
low-impedance line preprocessing yields smaller differences: 
0.0059 per unit voltage magnitude and 1.17° angle differ¬ 
ence. Thus, the differences from the preprocessing for the 
2869-bus PEGASE system can be largely attributed to the 
sensitivity of the OPE problem itself to small changes in the 
low-impedance line parameters. Preprocessing reduced the 
number of buses by between 21% and 27% for the PEGASE 
and between 9% and 26% for the Polish systems. 

These results show the need for further study of the sensi¬ 
tivity of OPE problems to low-impedance line parameters as 
well as additional numerical improvements of the moment 
relaxations and SDP solvers to reduce thrshz. 

V. Moment Relaxations and Penalization 

As will be shown in Section |Vl] the moment relaxations 
globally solve many large OPE problems with active power 
loss minimization objectives after removing low-impedance 
lines as described in Section m Directly applying the 
moment relaxations to many large OPE problems with more 
general cost functions has so far been computationally in¬ 
tractable. This section describes the non-convexity associated 
with more general cost functions and proposes a method 
to obtain feasible solutions that are at least near the global 
optimum, if not, in fact, globally optimal for many problems. 

Specifically, we propose augmenting the moment re¬ 
laxations with a penalization in the objective function. 
Previous literature [16], [17] adds terms to the first- 
order moment relaxation that penalize the total reac¬ 
tive po wer injection and the apparent p ower line loss 

(i-e., ifpim + fpmif + {fQini + fq-mif) for “problem¬ 
atic” lines identified by a heuristic. This penalization often 
finds feasible points that are at least nearly globally optimal. 

However, the penalization in [17] requires choosing two 
penalty parameters and fails to yield a feasible solution to 
some problems. This section describes progress in addressing 
these limitations by augmenting the moment relaxations 
with a reactive power penalization. The proposed “mo- 
mentH-penalization” approach only requires a single penalty 
parameter and finds feasible points that are at least nearly 
globally optimal for a broader class of OPE problems. This 
section also analyzes the convexity properties and provides 
a physical intuition for reactive power penalization. 

A. Penalization of Reactive Power Generation 

The penalization method proposed in [17] perturbs the 
objective function ([T]i to include terms that minimize the total 
reactive power loss and the apparent power loss on specific 
lines determined by a heuristic method. These terms enter 
the objective function with two scalar parameter coefficients. 
Obtaining a feasible point near the global solution requires 
appropriate choice of these parameters. 

For typical operating conditions, reactive power is strongly 
associated with voltage magnitude. Penalizing reactive power 
injections tends to reduce voltage magnitudes, which also 


tends to increase active power losses since a larger current 
flow, with higher associated loss, is required to deliver a 
given quantity of power at a lower voltage magnitude. 

For many problems for which the first-order moment 
relaxation fails to yield the global optimum, the relaxation 
“artificially” increases the voltage magnitudes to reduce 
active power losses. This results in voltage magnitudes and 
power injections that are feasible for the relaxation (ITTi) but 
infeasible for the OPF problem Q. 

By choosing a reactive power penalty parameter that bal¬ 
ances these competing tendencies (increasing voltage mag¬ 
nitudes to reduce active power losses vs. decreasing voltage 
magnitudes to reduce the penalty), the penalized relaxation 
finds a feasible solution to many OPF problems. Since losses 
typically account for a small percentage of active power 
generation and active and reactive power are typically loosely 
coupled, the reactive power penalization often results in a 
feasible point that is near the global optimum. 

We next study the convexity properties of the cost function 
and the reactive power penalization. The cost function ([T]) 
is convex in terms of active power generation but not 
necessarily in terms of the voltage components Vd and Vq due 
to the non-convexity of (iTaliPI Thus, the objective function is 
a potential source of non-convexity which may result in the 
relaxation’s failure to globally solve the OPF problem. 

Consider the eigenvalues of the symmetric matrices C 
and D, where, for the vector x containing the voltage 
components, is a linear cost of active power generation 

(i.e., C 2 = Co = 0 in ([T])) and x'^'Dx calculates the reactive 
power losses. For the 2383-bus Polish system [21], which has 
linear generation costs, the most negative eigenvalue of C 
is —8.53 X 10^. Thus, the objective function of the original 
OPF problem (fS^ is non-convex in terms of the voltage 
components, which can cause the relaxation (dB to fail 
to yield the global optimum. Conversely, active power loss 
minimization is convex in terms of the voltage components 
due to the absence of negative resistances. 

As indicated by the potential for negative eigenvalues of 
D (e.g., the matrix D for the 2383-bus Polish system has a 
pair of negative eigenvalues at —0.0175), penalizing reactive 
power losses is generally non-convex due to capacitive 
network elements (i.e., increasing voltage magnitudes may 
decrease the reactive power loss). See [18] for related work 
that uses a convex objective based on a Laplacian matrix. 

Further work is needed to investigate the effects of reactive 
power penalization on OPF problems with more realistic 
generator models that explicitly consider the trade-off be¬ 
tween active and reactive power outputs (i.e., generator “D- 
curves”). A tighter coupling between active and reactive 
power generation may cause the reactive power penalization 
to yield solutions that are far from the global optimum. 

The apparent power line loss penalty’s effects are not as 
easy to interpret as the reactive power penalty. Ongoing work 
includes understanding the effects of the line loss penalty. 

^The cost function of the moment relaxation GB is always convex. 
This section studies the convexity of the penalized objective function for 
the original non-convex OPF problem 



B. Moment+Penalization Approach 

Although the reactive power penalization often yields a 
near rank-one solution (i.e., rank-one matrices for most of the 
submatrices in the decomposition discussed in Section ITlI-Bl i. 
this penalization alone is not sufficient to obtain a feasible 
point for many problems. Reference [17] penalizes the ap¬ 
parent power line loss associated with certain lines to address 
the few remaining non-rank-one “problematic” submatrices. 
However, this approach has several disadvantages. 

First, penalizing apparent power line losses introduces 
another parameterO Introducing parameters is problematic, 
especially when lacking an intuition for appropriate values. 

Second, the combination of reactive power and line loss 
penalization may not yield a feasible solution to some 
problems. For instance, the OFF problems case9mod and 
case39modl from [1] are globally solved with low-order 
moment relaxations, but there is no known penalization of 
reactive power and/or apparent power line loss that yields a 
feasible solution for these problems. Also, the penalization 
approach is not guaranteed to yield a feasible solution that 
is close to the global optimum. 

Unlike the penalization approach, the moment relaxation 
approach does not require the choice of penalty parameters, 
globally solves a broader class of OFF problems, and is 
guaranteed to yield the global optimum when the rank 
condition ((TtTi is satisfied. However, direct application of the 
moment relaxations to large problems has so far been limited 
to active power loss minimization objective functions. We 
conjecture that the non-convexity associated with more gen¬ 
eral cost functions requires higher-order moment constraints 
at too many buses for computational tractability. 

To apply the moment relaxations to large OFF problems 
with active power generation cost objective functions, we 
augment the moment relaxations with a reactive power 
penalty. Specifically, we apply the sparsity-exploiting mo¬ 
ment relaxations described in Section [III] to the OFF prob¬ 
lem (|5]l where the objective function (l5al l is replaced by 

J2ifckiVd,Vg) + ebfQkiVd,Vg)) (13) 

k&g 

where ej, is the scalar reactive power penalization parameter. 
That is, rather than apply an apparent power loss penalization 
to the objective function, we apply higher-order moment 
constraints to specific buses [13]. As will be demonstrated in 
Section [Vl] higher-order moment constraints are only needed 
at a few buses in typical OFF problems after augmenting the 
objective function with a reactive power penalization term. 

Similar to the existing penalization, when the rank con¬ 
dition (IT^ is satisfied, the proposed “moment-tpenalization” 
approach yields the global solution to the modified OFF prob¬ 
lem but not necessarily to the original OFF problem (|5]). 
However, since the penalization does not change the con¬ 
straint equations, the solution to the moment-tpenalization 
approach is feasible for the original OFF problem ([T3I) . 

Reference [17] uses the same penalization parameter for each “problem¬ 
atic” line. Generally, each line could have a different penalty parameter. 


The first-order moment relaxation without penalization (i.e., 
= 0) gives a lower bound on the globally optimal 
objective value for the original OFF problem (|5]l. This lower 
bound provides an optimality metric for the feasible solution 
obtained from the momentH-penalization approach. As will be 
shown in Section [Vll the feasible solutions for a variety of 
problems are within at least 1% of the global optimum. 

The momentH-penalization approach inherits a mix of the 
advantages and disadvantages of the moment relaxation 
and penalization methods. First, the momentH-penalization 
approach requires selection of a single scalar parameter (one 
more than needed for the moment relaxations, but one less 
than generally needed for the penalization in [17]). This 
parameter must be large enough to result in a near rank- 
one solution (i.e., the solution should have many rank-one 
submatrices in the decomposition described in Section llll-BIl . 
but small enough to avoid large changes to the OFF problem. 

Second, the penalization eliminates the moment relax¬ 
ations’ guarantees: the momentH-penalization approach may 
yield a feasible solution that is far from the global optimum 
or not give any solution. However, the momentH-penalization 
approach finds global or near-global solutions to a broader 
class of small OFF problems than penalization approach 
of [17] (e.g., case9mod and case39modl with Cb — 0, and 
case39mod3 with Cb = $0.10/MVAr [1]). This suggests that 
the momentH-penalization approach inherits the ability of the 
moment relaxations to solve a broad class of OFF problems. 

Finally, the penalization in the momentH-penalization ap¬ 
proach enables calculation of feasible solutions that are at 
least nearly globally optimal for a variety of large OFF 
problems with objective functions that minimize active power 
generation cost rather than just active power losses. 

Note that it is not straightforward to compare the com¬ 
putational costs of the momentH-penalization approach and 
the penalization approach in [17]. A single solution of a 
penalized first-order moment relaxation, as in [17], is faster 
than a relaxation with higher-order moment constraints. 
Thus, if one knows appropriate penalty parameters, the 
method in [17] is faster. Although a relatively wide range 
of penalty parameters tends to work well for typical OFF 
problems, there are problems for which no known penalty 
parameters yield feasible solutions. For these problems, the 
momentH-penalization approach has a clear advantage. 

The momentH-penalization approach has the advantage of 
systematically tightening the relaxation rather than requiring 
the choice of penalty parameters. However, the higher-order 
constraints can significantly increase solver times. Thus, 
there is a potential trade-off between finding appropriate 
penalization parameters for the approach in [17] and in¬ 
creased solver time from the momentH-penalization approach. 
The speed of the momentH-penalization approach may be 
improved using the mixed SDF/SOCF relaxation from [14]. 

VI. Results 

This section first globally solves several large, active- 
power-loss minimizing OFF problems using moment re¬ 
laxations without penalization (cf, = 0). Next, this section 






applies the moment+penalization approach to find feasible 
points that are at least nearly globally optimal for several test 
cases which minimize active power generation cost. Unless 
otherwise stated, the preprocessing method from Section m 
with thrshz set to 1 x 10“^ per unit is applied to all 
examples. No example enforces a minimum line resistance. 

The results are generated using the iterative algorithm 
from [13] which selectively applies the higher-order moment 
relaxation constraints as summarized in Section IIII-BI The 
algorithm terminates when all power injection mismatches 
are less than 1 MVA. 

The implementation uses MATLAB 2013 a, YALMIP 
2015.06.26 [23], and Mosek 7.1.0.28, and was solved using a 
computer with a quad-core 2.70 GHz processor and 16 GB 
of RAM. The test cases are the Polish system models in 
Matpower [21] and several PEGASE systems [22] repre¬ 
senting portions of the European power system. 

A. Active Power Loss Minimization Results 

Table U shows the results of applying the moment relax¬ 
ations to several large OPE problems that minimize active 
power losses (i.e., the cost coefficients in ([T]i are C 2 = cq = 0, 
Cl = $l/MWh). The solutions to the preprocessed problems 
are guaranteed to be globally optimal since there is no 
penalization. The columns in Table U list the case name, 
the number of iterations of the algorithm from [13], the 
maximum power injection mismatch, the globally optimal 
objective value, and the solver time summed over all itera- 
tionsU The abbreviation “PL” stands for “Poland”. Table U 
excludes several cases (the 89-bus PEGASE system and the 
Polish 2736sp, 2737sop, and 2746wp systems) which only 
require the first-order relaxation and thus do not illustrate the 
capabilities of the higher-order relaxations. 


TABLE I 

Active Power Loss Minimization Results 


Case 

Name 

Num. 

Iter. 

Global Obj. 
Val. ($/br) 

Max 

(MVA) 

Solver 
Time (sec) 

PL-2383wp 

3 

24990 

0.25 

583 

PL-2746wop 

2 

19210 

0.39 

2662 

PL-3012wp 

5 

27642 

1.00 

319 

PL-3120sp 

7 

21512 

0.77 

387 

PEGASE-1354 

5 

74043 

0.85 

407 

PEGASE-2869 

6 

133944 

0.63 

921 


Each iteration of the algorithm in [13] after the first 
adds second-order constraints at two buses. Thus, a small 
number of second-order buses (between 0.1% and 0.7% of 
the number of buses in the systems in Table H] after the low- 
impedance line preprocessing) are applied to all examples 
in Table U This results in computational tractability for the 
moment relaxations. 

Note that PL-2746wop has a much greater solver time 
than the other systems even though it only has second-order 
constraints at two buses. This slow solution time is due to 
the fact that the two second-order buses are contained in 

^PEGASE-1354 and PEGASE-2869 use a thrshz of 3 X 10“® per unit. 
All other systems use 1 X 10“® per unit. 


Eigenvalue Ratio for PL-2383wp 



Fig. 1. Eigenvalue ratio (largest/second-largest eigenvalue) for each subma¬ 
trix for PL-2383wp. Large values (>10“^) indicate satisfaction of the rank 
condition EJ. For = 0 (green), most of the submatrices are not rank one. 
For = $50/MVAr (red), most submatrices satisfy the rank condition 
with the exception of those on the far left of the figure. Applying second- 
order moment constraints to two of the buses that are in these submatrices 
(blue) results in all submatrices satisfying the rank condition. 

submatrices corresponding to cliques with 10 and 11 buses. 
The second-order constraints for these large submatrices 
dominate the solver time. The mixed SDP/SOCP relaxation 
in [14] may be particularly useful beneficial for such cases. 

Since the low-impedance line preprocessing has been 
applied to these systems, the solutions do not exactly match 
the original OPE problems. Matpower [21] solutions of 
the original problems have objective values that are slightly 
larger than the values in Table Udue to losses associated with 
the line resistances removed by the preprocessing. 

After the low-impedance line preprocessing, local solu¬ 
tions from Matpower match the solutions from the moment 
relaxations and are therefore, in fact, globally optimal. This 
is not the case for all OPE problems [1], [13]. 

B. Moment+Penalization for More General Cost Functions 

As discussed in Section IV-AI minimization of active 
power generation cost often yields a non-convex objective 
function in terms of the voltage components. Despite this 
non-convexity, low-order moment relaxations typically yield 
global solutions to small problems, including problems with¬ 
out known penalty parameters for obtaining a feasible points 
(e.g., case9mod and case39modl from [1]). 

However, the moment relaxations are thus far intractable 
for some large OPE problems with non-convex objective 
functions. A reactive power penalty often results in the 
first-order moment relaxation yielding a solution that is 
nearly globally optimal (i.e., most of the submatrices in the 
decomposition described in Section IIII-BI satisfy the rank 
condition (ITSli i. Enforcing higher-order constraints at buses 
in the remaining submatrices yields a feasible solution to the 
OPE problem. This is illustrated in Eig. [T] which shows the 
ratio between the largest and second-largest eigenvalues of 
the submatrices of the moment matrix, arranged in increasing 
order, for the 2383-bus Polish system. If the submatrices 
were all rank one, then this eigenvalue ratio would be infinite. 
Thus, large numeric values (i.e., greater than 1 x 10"*) indicate 
satisfaction of the rank condition within numerical precision. 
Without the reactive power penalty, the rank condition is 
not satisfied for most submatrices. With the reactive power 
penalty, the rank condition is satisfied for many but not all 
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submatrices. Enforcing higher-order moment constraints at 
two buses which are in the high-rank submatrices results in 
a feasible (rank-one) operating point for the OPE problem 
which is within at least 0.74% of the global optimum. 

To further illustrate the effectiveness of the mo- 
mentH-penalization approach. Table HI] shows the results of 
applying the moment-npenalization approach to several large 
OPE problems with active power generation cost functions. 
The optimality gap column gives the percent difference 
between a lower bound on the optimal objective value from 
the first-order moment relaxation and the feasible solution 
obtained from the moment-tpenalization approach for the 
system after low-impedance line preprocessing. 

TABLE II 

Generation Cost Minimization Results 


Case 

Name 

($/MVAr) 

Num. 

Iter. 

Opt. 

Gap 

Max 

(MVA) 

Solver 
Time (see) 

PL-2383wp 

50 

2 

0.74% 

0.13 

152.2 

PL-3012wp 

50 

7 

0.49% 

0.20 

1056.3 

PL-3120sp 

100 

6 

0.92% 

0.08 

1164.4 


The penalized first-order relaxation requires 74.6, 88.9, 
and 97.0 seconds for PL-2383wp, PL-3012wp, and PL- 
3120sp, respectively. Attributing the rest of the solver time 
to the higher-order constraints implies that these constraints 
accounted for 3.1, 433.7, and 582.4 seconds beyond the time 
required to repeatedly solve the first-order relaxations. 

The momentH-penalization approach can yield feasible 
points that are at least nearly globally optimal for cases 
where both the penalization method of [17] and low-order 
moment relaxations fail individually. Eor instance, the mo¬ 
mentH-penalization approach with a reactive power penalty of 
€b = $0.10/MVAr gives a feasible point within 0.28% of the 
global optimum for case39mod3 from [1], but both second- 
and third-order moment relaxations and the penalization 
method in [17] fail to yield global solutions. 

VII. Conclusion 

“Moment” relaxations from the Lasserre hierarchy for 
polynomial optimization globally solve a broad class of OPE 
problems. By exploiting sparsity and selectively applying the 
computationally intensive higher-order moment relaxations, 
previous literature demonstrated the moment relaxations’ 
capability to globally solve moderate-size OPE problems. 
This paper presented a preprocessing method that removes 
low-impedance lines to improve the numerical conditioning 
of the moment relaxations. After applying the preprocessing 
method, the moment relaxations globally solve a variety of 
OPE problems that minimize active power losses for systems 
with thousands of buses. A proposed “momentH-penalization” 
method is capable of finding feasible points that are at least 
nearly globally optimal for large OPE problems with more 
general cost functions. This method has several advantages 
over previous penalization approaches, including requiring 
fewer parameter choices and solving a broader class of OPE 
problems. The results are demonstrated using large OPE 
problems representing European power systems. 
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